Strong Clustering of Faint Galaxies at Small Angular Scales 
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g ; ABSTRACT 

P_|. The 2-point angular correlation function of galaxies, uj{9) , has been computed on 

^ I equatorial fields observed with the CTIO 4m prime focus, within a total area of 2.31 

^ I deg^. In the magnitude range 19 < nin < 21.5, corresponding to < z >~ 0.35, we 

^ ■ find an excess of power in uj{9) at scales 2" < 6 < 6" over what would be expected 

from an extrapolation of uj{9) measured at larger 6. The significance of this excess is 
~ 5a. At larger scales, 6" < ^ < 24", the amplitude of oj{9) is 1.6 times smaller than 
^ ' the standard no evolutionary model. At these scales there is remarkable agreement 

■ between the present data and Infante & Pritchet (1995). 

At large angular scales (6" < 9 < 24") the data is best described by a model where 
clustering evolution in ^(r, z) has taken place. Strong luminosity evolution cannot 
be ruled out with the present data. At smaller scales, 2" < 9 < 6", our data are 
formally fit by models where e = — 2.4(f2 = 0.2, = 5.1/i~^Mpc) or = 7.3/i~^Mpc 
[Vt = 0.2, e = 0). If the mean redshift of our sample is 0.35 then our data show a clear 
detection of the scale (~ 19h~^kpc) where the clustering evolution approaches a highly 
non linear regime, i.e., e < 0. 

The rate at which galaxies merge has been computed. If this rate is proportional to 
{l + z)"", then m = 2.2 ± 0.5. 

Subject headings: galaxies: evolution - galaxies: format ion - galaxies:interactions 



^CNPq Fellow 



1 



1. Introduction 

The two point angular correlation function 
has been extensively used to study the clustering 
properties of galaxies. At large angular separa- 
tions, 9 > 10", it is well established that faint 
galaxies are less clustered compared to galaxies 
observed locally. Infante &; Pritchet (1995) using 
observed rcdshift distributions and straightfor- 
ward calculations with little model dependency 
found that predicted angular correlation ampli- 
tudes are about a factor of two higher than those 
observed over a wide range from 20 < < 24. 
The blue galaxy excess in the counts N(m), and 
the correlation deficiency, can be neatly resolved 
by introducing a "new", weakly clustered popu- 
lation at low redshifts of 0.2 < z < 1 (Efstathiou 
et al. 1991; Infante k Pritchet 1995). However, 
luminosity-dependent evolution must be invoked 
in order to brighten up invisible low-L galaxies 
at z = to ~ L* at z ~ 0.3, while still maintain- 
ing the no-evolution shape of N(z) (Colless et al. 
1990; Lilly et al. 1991; Broadhurst et al. 1988). 

At small angular scales, 6 < 10", the situa- 
tion is different and there is no consensus yet. 
A limited number of studies have been success- 
ful in measuring the angular correlation function 
at small separations {9 < 6", physical separa- 
tion less than ~ 20h~^kpc at z ~ 0.3)(Carlberg 
et al. 1994; Ncuschaefcr et al. 1995; Woods et 
al. 1995). The main difficulty in such determi- 
nations is poor statistics on scales < 10". The 
results range from a strong excess to a deficiency 
of pairs with respect to an inward extrapolation 
of io{9) oc 9~^'^ at larger separations. This excess 
in u!{9) may be due to a population of "compan- 
ions" not present at the current epoch, or lumi- 
nosity enhancement of intrinsically faint galaxies 
in pairs. 

In this letter we present the first clear detec- 
tion of a break from a power law observed at 
larger separations in the angular two point cor- 
relation function, uj{9), at 9 < 6". We have mea- 
sured uj{9) at separations 2" < < 68" on a 2.31 



deg2 field at 19 < tur < 21.5 (< z >^ 0.35) from 
a recently selected catalog of faint galaxy pairs 
and groups (de Mello et al. 1996a and b, here- 
after Paper I and III). Special care was taken in 
this work to resolve galaxies at small separations. 

2. The Data 

The data set is comprised of 96 equatorial im- 
ages of 15' X 15' (0.44" pixel"^) making a total 
area of 5.25 deg^ taken at the CTIO 4m prime 
focus camera by the Higli-Z Supernovae Search 
Group (Lcibundgut ct al. 1995) over two nights 
in March and November 1995. The fields are all 
equatorial within O'' < R.A. < 4:^ and lO'' < 
R.A. < 14^. Center field coordinates will be 
given in a forthcoming paper. Five minute ex- 
posures through a redshifted B filter which is al- 
most equivalent to a regular Kron- Cousin R filter 
(B/(z=0.4), hereafter mn) were sufficient to pro- 
vide good quality images. After a careful inspec- 
tion of the images we decided to run our analysis 
on the best images. Our selection criteria were 
based on seeing, low number of bright stars and 
photometric conditions. The total eff'ective angu- 
lar area on which we have computed the angular 
correlation function is 2.31 deg^. 

The observations, reductions, calibrations and 
selection criteria are described in Paper I. Coor- 
dinates of each pair /group is also given in Pa- 
per I and III. A number of tests were conducted 
in order to determine detection performance and 
limiting magnitudes. We claim in Paper I that 
our catalog is 99% complete at m/j < 22, 
and that the two components of a pair sepa- 
rated by more than 2" were always detected for 
1.0" < seeing < 1.6". Star/galaxy separation 
was done using Kron photometry (Kron 1980) 
and the properties of the inverse first and second 
moments of the images which gives a measure of 
intrinsic size and central compactness (see Paper 
III for more details). Results were tested by eye 
inspection. 
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3. Analysis and Results 

The angular correlation function, uj{6), was 
computed using a method that is more fully de- 
scribed in Infante (1994; see also Infante and 
Pritchet 1995, hereafter IP95). An artificial cat- 
alog of randomly distributed objects was created, 
and areas contaminated by defects, bright stars, 
etc. were masked from both the artificial catalog 
and the real catalog. The estimator used was 



NggNr 



B Ngr {Ng - 1) Nrr, {Nr " 1 



where Ngg is the number of galaxy pairs in a 
given range of separations summed over all CCD 
fields, Ngr is the number of random pairs using 
galaxies as centers, Ng is the number of galaxies, 
A^,. is the number of random objects, Nrr is the 
number of random pairs, and Nm is the same 
as Nrr but an independent set of random counts. 
B is the "integral constraint" correction factor 
(e.g. Peebles 1980; Koo and Szalay 1984). At 
the relevant scales in this paper, 6 < 24", the 
uncertainties are limited by Poisson noise rather 
than by variations in B. We have used B = 1 for 
our clustering at small separation analysis. 

The angular correlation function uj{9) was de- 
rived in a 2.31deg^ area. Fig. || gives the an- 
gular correlation function at separations 9 < 2' , 
for the magnitude interval 19< mji <21.5. Also, 
the canonical lo=2.39~^'^ , as derived in IP95, 
is given for comparison in Fig. |^ (thick solid 
line). The uncertainties, 6oj, are 65% confidence 
intervals computed using the bootstrap resam- 
pling method as described in Efron &; Tibshirani 
(1986). Five hundred resamplings were carried 
out in each computation. We note that Poisson 
errors are ~ 2.5 times smaller than bootstrap er- 
rors. In most data points the 2a Poisson error 
bars are smaller than the symbol. 

Using least squares to fit uj{9) we obtain u = 
(2.06 ± 0.02)6l-° '^^^° °2, in the range 6" < 9 < 
67.6". A K-S significance test shows that in this 



range, both the canonical IP95 and lo{9) for our 
data, are drawn from the same distribution with 
a 90% significance level. 

Table || shows the number, N, of pairs at sep- 
arations 2" < 9 < 6" and 6" < 9 < 24" for the 
real and random catalogs. The fractional excess. 
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in the number of pairs with respect to the random 
catalog are 0.637 and 0.200, respectively, while 
(4^) over what would be expected from 

\ " / pairs 

IP95 are 0.361 and -0.011, respectively. We also 
show in Table || the number of isolated groups, 
with 2, 3, 4, 5 and 6 members, found in our cat- 
alog and the number of groups found in the ran- 
dom catalog. 

In Fig. ^ the fractional excess number of pairs, 
(^) , for 19 <mR< 21.5 are plotted with 

V / pairs 

respect to the random catalog expectation (thin 
line), to IP95 data (thick line) and to various 
models (dashed lines) as explained in the next 
section. At scales 9 < 6" our data departs from 
the canonical uj{9) from IP95 by ~ 5a. This ex- 
cess power has not been detected in most previ- 
ous works, only Carlberg et al. (1994) report an 
excess in the number of pairs at separations < 6" 
and V < 22.5. 

4. Discussion 

Are there any effects in the data or data re- 
duction techniques that might in principle cause 
an excess in uj{9) at small separations? First, 
we consider the possibility that rich clusters of 
galaxies at z ~ 0.3 in our fields might signif- 
icantly enhance uj{9) at 2" < 9 < 6". From 
Table || the contribution of groups with more 
than 5 members is 52 pairs corresponding to a 
negligible (^) = 0.032. Next, we con- 

V ^' / pairs 

sider multiple detections of several faint spurious 
sources for each physically distinct galaxy by our 
finding algorithm. For this purpose, all images 
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were inspected by eye, eliminating all possible 
sources of contamination, namely, bright galax- 
ies, scratches, bad pixels and cosmetic defects on 
the chip. There is also the possibility of contam- 
inating light. Galaxies that normally would have 
been below the magnitude limit may be raised 
above it by the light of the second object. Our 
photometry accounts for this effect. We do an 
excision of disturbing objects for crowded field 
photometry. Our algorithm takes the median in 
rings and replaces deviant values with the me- 
dian. And finally, we checked for systematic 
effects on magnitude determination. Extensive 
simulations were performed in order to test the 
detection and photometry of faint images as a 
function of magnitude. No systematic effects 
were detected. See Paper I for details regarding 
the simulations. 

4.1. Comparison with other Observations 

Several studies have investigated io{9) of faint 
galaxies in a number of passbands (IP95 and ref- 
erences therein; Neuschaefer et al. 1995; Shep- 
herd et al. 1996). The principal result from 
these surveys is that uj{9) follows a power law, 
uj = A e-\ where 0.7 < 6 < 0.9 {6 = 0.8 
is the standard value), and A is the amplitude 
scaling with magnitude. In this work we have 
taken IP95 as the canonical io{9) in the red and 
the solid line in Fig. |l] represents this power 
law, Lo = 2.3 9~^'^. After transforming IP95 
Fband to m^ of this work (F=mij-F0.094(B-R), 
(B-R)=1.68, Metcalfe et al. 1991) we find that 
the agreement for 6" < ^ < 68" is remarkable. 
However, at small angular scales {9 < 6"), except 
for the work by Carlberg et al. (1994), no pre- 
vious work shows a significant excess in io{9) by 
more than one a. The quality and amount of 
data that we have used in this work are sufficient 
to make us believe that the excess at small scales 
is a real effect. 



4.2. Comparison with Models 

How can a high clustering amplitude at small 
separations be understood? One of the possibil- 
ities is to assume that mergers were more fre- 
quent at z « 0.35 and the excess of pairs that 
we are seeing at that epoch is in fact mergers 
in progress. They represent a true density ex- 
cess. Another scenario that we cannot rule out 
is the possibility of an excess of strongly clustered 
dwarf galaxies. Two alternatives are possible in 
this hypothesis: a new class of dwarf population 
proposed by Babul & Rees (1992) which would 
be recently 'active', but subsequently fading out 
of view; or the 'sub-units' proposed by Broad- 
hurst et al. (1992) which would slowly merge to 
form more massive galaxies. It is also possible 
that this excess in uj{9) is due to brightening of 
low luminosity companions (Carlberg, 1992). A 
fourth possibility is that although there is clus- 
tering evolution at scales larger than 6" there is 
no clustering evolution at scales < 6". In this 
case, CO (9) is well fitted by a model with e = at 
2" <9 < 6" and e = 0.8 at 6* > 6". Thus, there is 
a scale 9 « 6" (~ 20h~^kpc, z ~ 0.35) at which 
the clustering evolution slows down. 

In order to show this latter point we have com- 
puted standard models for uj{9) . If we assume a 
power law spatial correlation function, 

e(r,z)= (^^) \l+z)-(3+^), (3) 

where r is the proper distance, is the proper 
correlation length and e is the clustering evolu- 
tion index; e = 7 — 3 means clustering fixed in 
comoving coordinates and e = represents stable 
clustering in physical coordinates. The angular 
correlation function is obtained through Limber's 
equation (Peebles 1980 eqs [56.7] and [56.13] and 
IP95). The no-evolution redshift distribution, 
is taken from Metcalfe et al. (1991) (see IP95 
for discussion and justification of our choice). 

In Fig. H three models are shown: The stan- 
dard model (7 = 1.8, e = 0, = 0.2 and 
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To = 5.1 h~^Mpc, solid line), a model with clus- 
tering evolution, (7 = 1.8, e = 0.8, Q = 0.2 and 
To = 5.1 h~^Mpc, dashed line), and a model with 
either clustering evolution, (7 = 1.8, e = —2.4, 
Q = 0.2 and To = 5.1 h~^Mpc, dashed dot- 
ted line) or a model with a strong local corre- 
lation length (7 = 1.8, e = 0, $7 = 0.2 and 
To = 7.3 h~^Mpc, dash-dotted line). 

The standard model with clustering evolution 
(e = 0.8) fits very well the observations in the 
range 6" < 9 < 24". The value of e agrees with 
the results obtained by Shepherd et al. (1996), 
where the low amplitude of the correlation func- 
tion might be due to strong clustering evolution 
since z ^ 0.35. 

The small scale data {6 < 6") is better fit- 
ted by a model with either a high local cluster- 
ing length (ro = 7.3 h'^Mpc) or e = -2.4. If 
we assume a correlation length as observed for 
present day L* galaxies (ro = 5.1 h~^Mpc), no 
clustering evolution (e = 0), and no luminosity 
evolution, the clustering properties of galaxies at 
< z >= 0.35 and spatial separations ~ 20h~^kpc 
(for z 0.35), in the niR passband are similar to 
present day L* galaxies. The excess in u{6) with 
respect to this model might correspond to the 
galaxy merger rate discussed in the next section. 

If the mean redshift of our sample is 0.35 then 
our data show a clear detection of the scale, 

19h~^kpc, at which the clustering evolution 
approaches a highly non linear regime, where 
e < is expected. 

4.3. The galaxy merger rate 

The merger rate is expected to have a strong 
redshift dependence varying as roughly as (1 -|- 
z)"^, where m ~ 4.5r2^'^. However, the diffi- 
culty of identifying objects even at moderate red- 
shift is a limiting factor in studies of faint galaxy 
pairs. Another factor is imposed by the signifi- 
cant amount of large telescope time required in 
order to obtain redshifts and good quality im- 
ages for a statistically significant sample of faint 



galaxies. So far, only a few previous works have 
redshift information for faint galaxies in pairs 
(see Koo k Kron 1992 and Elhs 1995 for re- 
views on redshift surveys and imaging of the faint 
galaxy population). 

There are two ways (not independent, though) 
to calculate the galaxy merger rate from our 
data. First, we count the number of isolated 
pairs of galaxies in the magnitude range 19 < 
mR < 21.5, corresponding to < 2; >^ 0.35, at 
scales 2" < 9 < 6" {^Ah-^kpc - 19.3hr^kpc, for 
z 0.35 and qo = 0.1), and compare it with 
numbers of pairs from local galaxy surveys se- 
lected in the same way. The fraction of galaxies 
in pairs in our catalog is 10.04ib0.17% {^/n uncer- 
tainty). The expected fraction of physical pairs is 
0.636 ±0.007, if a;(2"-6") = 1.75 ±0.05. Follow- 
ing Carlberg et al. (1994), the fraction of galax- 
ies that will merge is 0.7. Thus, the percentage 
of galaxies in physical pairs that will eventually 
merge is 10.04 x 0.636 x 0.7 = 4.47 ± 0.01%. At 
low redshifts this value is 2.3 ± 0.2% (Carlberg 
et al. 1994). If the rate at which galaxies merge 
is proportional to (1 + z)*", then m = 2.2 ± 0.5. 
(The uncertainty in the merger rate is dominated 
by the uncertainty in the fraction of galaxies that 
will merge. We assume 0.1.) 

Second, the galaxy merger rate can be cal- 
culated using the excess correlation of our faint 
galaxies over what would be expected from the 
standard model. The fractional excess ( ^ff- ) 

V /pairs 

over what would be expected from the standard 
model is 0.332 (< 9 >= 4" corresponding to 
15.1h~^kpc at < z >= 0.35). Given that 
(1 + z)"" = ^Yien m = 1.34. It is 

^ ' _ i-+<^standard 

to be noticed that no corrections for the fraction 
of physical pairs that will merge has been made 
in this case. 

5. Conclusions 

We have estimated uj{9) from a well defined 

sample of galaxies at 19 < niR < 21.5 and at 
small angular separation, 2" < 9 < 6". We have 
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performed simulations in order to test detection 
efficiency of galaxy pairs and the limiting mag- 
nitude in our galaxy sample. We find that our 
algorithms can detect, at a 99% confidence level, 
the two components of a pair separated by > 2" 
and with < 22.5 in our CTIO 4m prime focus 
CCD images. 

At &' < 6 < 68" the agreement between the 
present data and IP95 is remarkable. The am- 
plitude of uj{9) at these separations is 1.6 times 
smaller than the no-evolution standard models. 

Our analysis indicates that there is an excess 
of power in oj{6) at 2" < 9 < 6" over what would 
be expected from an extrapolation of the canon- 
ical power law uj{9) at larger 9. The significance 
of this excess is ~ 5a. These results place limits 
on the galaxy merger rate of (l+z)^-^~^-^. 

At large angular scales (6" < 9 < 24") the 
data is best described by a model where cluster- 
ing evolution in ^(r, z) has taken place. Strong 
luminosity evolution cannot be ruled out with 
the present data. At smaller scales, 2" < 9 < 
6", our data are formally fit by models where 
e = —2.4(0 = 0.2) or the clustering length 
ro = 7.3h-^Mpc. If the mean redshift of our 
sample is 0.35 then our data show a clear detec- 
tion of the scale (~ 19h~^kpc) where the clus- 
tering evolution approaches a highly non linear 
regime, i.e., e < 0. 

We are grateful to the High-Z Supernovae 
Search Group for making their images avail- 
able. DFM thanks CNPq for the fellowship and 
LI thanks Fondecyt Chile for support through 
'Proyecto 1960414'. 

REFERENCES 

Babul, A., & Rees, M. J. 1992, MNRAS, 255, 346 

Broadhurst, T. J., Elhs, R. S. & Glazebrook, K. 
1992, Nature, 355, 55 

Broadhurst, T. J., Ellis, R. S., & Shanks, T. 1988, 
MNRAS, 235, 827 



Carlberg, R. G. 1992, ApJ, 399, L31 

Carlberg, R. G., Pritchet, C. J., & Infante,L. 
1994, ApJ, 435, 540 

Colless, M., Ellis, R. S., Taylor, K., & Hook, R. 
N. 1990, MNRAS, 244, 408. 

de Mello, D. F., Infante, L., &; Menanteau, F. 
1996a, ApJSS submitted, Paperl 

de Mello, D. F., Infante, L., &; Menanteau, F. 
1996b, in preparation, Paperlll 

Efron, B., & Tibshirani, R. 1986, Statistical Sci- 
ence, 1, 1, 54 

Elhs, R. S. 1995, in Stellar Populations, lAU 
Symposium 164, eds P. C. van der Kruit &i 
G. Gilmore (Kluwer), p. 291 54 

Efstathiou, G., Bernstein, G., Katz, N., Tyson, 
J. A., k Guhathakurta, P. 1991, ApJ, 380, 
L47 

Infante, L. 1994, A&A, 282, 353 

Infante, L., & Pritchet, C. J. 1995, ApJ, 439, 565. 
IP95 

Koo, D. C, & Kron, R. G. 1992, ARA&A, 30, 
613 

Koo, D. C, & Szalay, A. S. 1984, ApJ, 282, 390 

Kron, R. 1980, ApJSS, 43, 305 

Leibundgut, B., et al. 1995, The Messenger, 81, 
19 

Lilly, S. J., Cowie, L. L., k Gardner, J. P. 1991, 
ApJ 369, 79 

Metcalfe, N., Shanks., T., Fong, R., & Jones, L. 
R. 1991, MNRAS, 236, 207 

Neuschaefer, L. W., Ratnatunga, K. U., Griffiths, 
R. E., Casertano, S., & Im, M. 1995, ApJ, 453, 
559 



6 



Peebles, P. J. E. 1980, The Large-Scale Structure 
of the Universe. Princeton University Press. 

Shepherd, C. W., Carlberg, R. G., Yee, H. K. C, 
k, Elhngson, E. 1996, preprint 

Woods, D., Fahlman, G. G., & Richer, H. B. 
1995, ApJ, 454, 32 



This 2-columii preprint was prepared with the AAS 
macros v4.0. 



7 



Table 1 
Excess Number of pairs. 



Catalog 


Oi-92 
[arcsec] 




^total 


j^obs 
pairs 


j\jran 
pairs 


2 


3 


Groups 
4 


5 


6 


(^) . ' 

\ ' pairs 


Real 


2-6 


1.753 


16749 


1317 


488 


841 


114 


15 


4 


3 


0.361 


Random 


2-6 


-0.010 


16749 


477 


497 


438 


14 













Real 


6-24 


0.250 


16749 


9962 


7921 


1318 


314 


58 


13 


2 


-0.011 


Random 


6-24 


0.004 


16749 


7970 


7943 


1587 


255 


30 


5 


1 





^Fractional excess number of pairs over what would be expected from IP95. 
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Fig. 1. Angular Correlation Function for 
19 < rriR < 21.5. The open circles are the data. 
65% uncertainties (500 bootstrap resampling of 
the data) are shown as vertical lines. The dark 
sohd hue is the IP95 lj{6) result with 5 = —0.8. 
The thin solid line is the no-evolution model de- 
scribed in the text (ro = 5.1 Mpc, e = 0). 
The dashed line is a model with clustering evo- 
lution (ro = 5.1 h-^ Mpc, e = 0.8). The dashed 
dotted line is a model with either clustering evo- 
lution (ro = 5.1 Mpc, e = —2.4) or strong 
correlation length (ro = 7.3 /t"^ Mpc, e = 0). 

Fig. 2. — Fractional excess number of pairs for 
19 < ruR < 21.5, as defined in equation 2. The 
dark solid line is the IP95 lo{6) result with S = 
—0.8. The dashed line is the no-evolution model 
described in the text (rg = 5.1 Mpc, e = 0). 
The thin solid line is the fractional excess with 
respect to a random distribution of points. 
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